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ABSTRACT 

We study the effects of rotation on the neutrino emission from Population 
III (Pop III) stars by performing a series of two-dimensional rotational collapse 
simulations of Pop III stellar cores. Our results show that rotation enhances 
the neutrino luminosities and the average energies of emitted neutrinos. This is 
because the thermalized inner core, which is the dominant neutrino source from 
Pop III stars, can be enlarged, due to rotational flattening, enough to extend 
the inner core outside the neutrinospheres. This is in sharp contrast to the case 
of spherical collapse, in which the case of inner core shrinks deeper inside the 
neutrinospheres before black hole formation, which hinders the efficient neutrino 
emission. In the case of rotational core-collapse, the emitted neutrino energies 
are found to become larger in the vicinity near the pole than the ones near the 
equatorial plane. These factors make the emergent neutrino spectrum broader 
and harder than the spherical collapse case. By computing the overall neutrino 
signals produced by the ensemble of individual rotating Pop III stars, we find that 
the amplitudes of the relic neutrinos, depending on their star formation rates, can 
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dominate over the contributions from ordinary core-collapse supernovae below a 
few MeV. A detection of this signal could be an important tool to probe star 
formation history in the early universe. 

Subject headings: hydrodynamics — methods: numerical — neutrinos — stars: 
rotation 



1. Introduction 



One of the most important goals in modern cosmology is to understand how the first 
stars formed at the end of the dark ages, and how they transformed the initially simple, ho- 



mogeneous universe into a state of ever increasing complexity (e.g., iBarkana &; Loebl 12001 



Bromm fc Larson! 12004 : lLoebl 120061 ) . The first stars, so called Population HI (Pop HI), are 



predicted to have been predominantly very massive with M > lOOM^ fe.g. JNakamura Sz Umemura 



2001: Abel et al. 2002; Bromm et al. 2002; Yoshida et al. 2006; O'Shea fc Norman 



2003). 



During their evolutions, the central core is thought to experience an electron-positron pair 
creation instability after carbon burning, which reduces the thermal energy of the core and 
eventually triggers gravitational core-collapse. For stellar masses less than ~ 26OM , rapid 
nuclear burning releases a large amount of energy sufficient to entirely disrupt the star as 
pair-instability supernovae. More massive stars, which also encounter pair-instability, are so 
tightly bound that the fusion of o xygen is unable to reverse infall. S uch stars are thought 



to collapse into black holes (BHs) (IBond et al.lll984l ; iFryer et al.ll200ll ). which we investigate 



in this paper. Information about the formati on and evolut i on of the Pop III stars might be 
observationa l ly ob tained from reionization ( Alvarez et al. 2006h . the infrared backgroun d 



gamma-ray bursts ( 


Schneider et a 


2002 


(Buonanno et al. 


20051; 


Suwa et al. 


2007a) 



yHgOod ; llwamoto et all 120051 ) 



20021 ; iBromm fc Loebl l2006h . and gravitational 



waves 



2007al ). However at present, it is still unclear whether or 



not one can obtain clear information about Pop III stars directly from observations. 

Recently, neutrinos, which should have been emitted in large numbers at the phase of 
core-collapse, are expected to be new eyes to unveil the mysteries of the first stars. Since 
the direct detection is unlikely because of their large distances, a diffuse background as 
relic neutrinos could potentially be observable. In an analogous context, the relic neu- 
trino background from ordinary c ore-collapse supernovae have been s t udied elaborately 



2003; 



(IBisnovatvi-Kogan fc Seidovl Il984l ; iTotani et al.l Il996l ; iKaplinghat et al.l |2000| ; lAndo et al. 



Ando fc Sato! 2004 ). Since the physics involved in relic neutrinos covers a very wide 
range of topics, e.g., cosmic star formation rate and the collapse dynamics, detecting the 
relic neutrinos or even setting limits on their flux are expected to give us quite useful and 
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unique implications not only for the understanding of the collapse physics itself, but also for 
various fields of astrophysics and cosmology. 

Th e detectabili t y of t he relic neutrinos from Pop III stars has been discussed by several 



papers. Ilocco et al.l (120051 ) estimated the amplitudes of the relic neutrinos produced during 
the nuclear burning phases as well as the core collapse phase, assuming a bar yonic fraction 
of Pop III stars of 10~ 3 and monochromatic progenitors with a mass of 300M Q . iDaigne et al. 



(120051 ) investigated the relic neutrino background produced by an early burst of Pop III stars 
taking into account a "normal mode" of the star formation at low redshift. The framework 
of hierarchical structure formation, on which their computation relies, is based on the cosmic 
star formation histories, which are constrained by the observed star formation rate at redshift 
z < 6, the observed chemical abundances in damped Lyman a absorbers and in the IGM, 
and allow for an early reionization of the Universe at z ~ 10 — 20. Regardless of the 
differences in the treatment of the star formation histories, these studies concluded that the 
typical energy of the relic neutrino background is lowered as low as MeV or sub-MeV due to 
the cosmological redshift so that the detection is out of the question with presently known 
experimental techniques. However it should be noted that in their studies, the neutrino flux 
and energy emitted from the Pop III stars are set by hand due to the lack of Pop III star 
collapse simulations. 

For more precise estimates of relic neutrinos from Pop III stars, we have to perform the 
hydrodynamical simulations. So far, there have been only a few simulations studying the 
gravitational co llapse of Pop HI stars with BH formation including appro pr iate microphys 



ical treatment ( Fryer et al. 



200ll ; iNakazato et all 12009 : ISuwa et al.ll2007bf ). iNakazato et al. 



(120061 ) performed a spherically symmetric general relativistic simulation in a wide range of 
masses (300 — 10000M o ), in which the state-of-the-art neutrino physics are taken into ac- 
count. Their detailed calculations revealed the properties of the emergent neutrino spectrum, 
and based on that, they estimated the relic neutrino background and found that the detection 
for the spherically collapsing Pop III stars is difficult for the currently operating detectors. 
Meanwhile current wisdom tells us that very metal poor stars lose very small amounts of 
mass and angular momentum through radiati yely driven stellar winds so that Pop III stars 
may collapse with large angular momentum (IHeger et al.l 120031 ) . So stellar rotation might 
vary the spectrum of emi tted neutrinos and affect its potential observability, which is the 
motivation of this study. iFryer et al.l (120011 ) investigated the collapse of a rotating Pop III 
star of mass 3OOM with gray neutrino transport simulations and discussed effects of ro- 
tation on the emitted luminosities. They, however, investigated only one rotational model. 
Considering the uncertainty of the angular momentum distributions of the Pop III stars, a 
more broad study is required. More recently, we implemented a ser ies of two-dimensi onal 
magnetorotational core-collapse simulations of Population III stars (ISuwa et al.ll2007bl ). in 
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which our interest was placed on the MHD effects on the formation of the jets. 

In this study, we focus on the effects of rotation on the neutrino emission during the 
core-collapse phase of Pop III stars. For this purpose, we carry out two-dimensional simula- 
tions of rotational core-collapse of Pop III stars, changing the initial strength of the angular 
momentum in a parametric manner. By so doing, we systematically investigate the depen- 
dence of the luminosity, the average energy, and the spectrum of neutrinos upon the rotation 
strength. We also pay attention to the anisotropy of the emergent neutrino energies, that 
is, the neutrino emission at polar direction and in the equatorial plane. Using the spectrum 
of single Pop III stars, we calculate the number flux of the relic neutrino background by 
summing up the contribution from individual Pop III stars and discuss their detectability. 

This paper is organized as follows: In the next section, we briefly introduce our numerical 
methods and the initial models. In §3 and §4, we present the numerical results. We start 
with the outline of the dynamics and neutrino emission in the case of spherical collapse 
(§3]) and then move on to discuss the effects of rotation on the dynamics and the neutrino 
luminosity ( §4.11) . the neutrino energies ( §4.2p and the spectrum ( §4.3p . §5 is devoted to 
summary and discussion. 



2. Numerical Techniques and Initial Models 



As for the hydro solver, we employ the ZEUS-2D code (IStone &: Norman! Il992l ) as a 
base and added major changes to include the microphysics appropriate for the simulation 
of Pop III stars. First we have added an equation for the electron fraction in order to 
treat electron and positro n captures and have approximated the neutrino transport by the 



so-ca l led leakage scheme ([Ruffert et al.l Il996l ; iKotake et al.l 120031 ; iRosswog fc Liebendorfer 



20031 ; iTakiwaki et al.l 120071 ). We consider three neutrino flavours: electron neutrinos, 



electron anti-neutrinos, u e , and the heavy- lepton neutrinos, u^, z? M , i/ T , v T , which are collec- 
tively referred to as fx- Reactions of u x , pair, photo, and plasma processes are included 
using the rates by lltoh et al.l (119891 ). As for the equation of state, we have incorporated the 
tabulated one based on relativistic mean field theory (see, IShen et al.lll998l ). In our 2D cal- 
culations, axial symmetry and reflection symmetry across the equatorial plane are assumed. 
Spherical coordinates (r, 9) are employed with logarithmic zoning in the radial direction and 
linear zoning in 9. One quadrant of the meridian section is covered with 300 (r) x 30 (9) mesh 
points. We also calculated some models with 60 angular mesh points, however, there was no 
significant difference. Therefore, we report only the results obtained from the models with 
30 angular mesh points. General relativistic gra yity is taken into ac count approximately by 
an "effective relativistic potential" according to iMarek et al.l (120061 ) . which is newly added 
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in our code. 

In this paper, we set the mass of the Pop III star to be 3OOM . This is consistent with the 
recent simulations of the star-formation phenomena in a metal free environment, providing 



Nakamura fc Umemura 



an in itial mass function peaked at masses 100 — 300M Q (see, e.g., 
200ll ) . We choose the value because we do not treat the nuclear-powered pair instability 
supernovae (M < 26OM ) and also t o facilitate the co mparison with the previous study, 
which employed the same stellar mass (IFryer et al.ll200ll ). 



We start the collapse simulations with a 180M Q core for the 300M© star. The core, 
which is the initial co ndition of ou r simu lations, is produced in the following way. According 
to the prescription in iBond et al.l (119841 ). we set the polytropic in dex of the core t o n = 3 
and assume that the core is isentropic with ~ 10ks per nucleon (IFryer et al.ll200ll ) with a 
constant electron fraction of Y e = 0.5. We adopt a central density of 5 x 10 6 g cm -3 at 
which the temperature of the central regions become high enough to photodisintegrate the 
iron (~ 5 x 10 9 K), thus initiating the collapse. Given the central density, the distribution 
of electron fraction, and the entropy, we numerically construct the hydrostatic structures of 
the core. 

Since we know little about the angular momentum distribution in the cores of Pop III 
stars, we assume the shellular rotation follows the initial rotation law: 



fi(r) = tt 



r 2 + r. 



2 • 



where Q(r) is an angular velocity, r is a radius, and Qo and r are model parameters, which 
determines a total angular momentum and the degree of differential rotation, respectively. 
In this paper, we fixed the value of r as 10 9 cm. Since the radius of the outer edge of the core 
is taken to be as large as 3.5 x 10 9 cm, the above profile represents a mildly differentially 
rotating core. Changing the initial rotational energies by varying the values of Qq, we 
compute 8 models, namely, one spherical and 7 rotational models. By so doing, we hope to 
see clearly how the dynamics and the neutrino emission deviate from those in the spherical 
case. We change the initial values of T/IW 7 ! from 0.01 to 2 %, where T/|W| represents the 
ratio of the rotational to the gravitational energy. We take the model with T/|W| 0.5 % 
as the canonical rotating model. This is because the rotation rate is tak en from the stellar 
evolution calculations of extremely rapidly rotating cores of massive stars (IHeger et al.ll2000i ) 
pushed by the study that Pop HI stars cou ld rotate rapidly due to insufficient mass-loss in 
the main sequence stage (IHeger et al.ll2003l ). 



The shape of the neutrino spheres is important for analyzing the properties of neutrino 
emission. Following the common practice, we define the neutrinosphere, R u , satisfies the 
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(2) 



condition 

" dr _ 2 
T~ 3' 

where A is the mean free path of neutrinos. We perform the integration for each angular bin 
to obtain R v {9). It should be noted that the notion of neutrinosphere is rather ambiguous. 
Neutrino reactions are highly energy dependent, and so the neutrinosphere should be as well. 
In theory, we should d istinguish the last scattering surface from the surface of last energy 
exchange (|Jankalll995l ). However, we stick to this approximate notion here for simplicity. 



3. Spherical Collapse 

We first outline collapse dynamics in the case of spherical symmetry. In Figure (U we 
show the time evolution of density (top left), temperature (top right), and radial velocity 
(bottom left) as a function of radius. Also shown is the radial velocity as a function of the 
mass coordinate (bottom right). We show the results for five time slices before BH formation 
(see figure for detail). As indicated in this figure, the core collapse of a very massive star can 
be divided into two phases; the infall phase (—24.1 and —10.1 ms) and the accretion phase 
(from —2.0 to ms). 

The infall phase sets in due to the gravitational instability induced by the photodisinte- 
gration of iron and the neutronisation that have the effect of reducing the pressure support 
at the core of the star, triggering a rapid collapse of the core just as an ordinary supernovae 
progenitor. However, there are several important differences between the structures of very 



massive stars and ordinary supernovae . As noted by previous works (e.g.. lBond et al.lll984 



Fryer et al.l |2001| ; iNakazato et al.l 120071 ) , the entropy in the more massive core is larger and 
thus favors more complete photodisintegration of heavy elements and a-particles so that 
the source of instability is reduced. This leads a thermally stabilized core □ unlike ordinary 
supernovae, whose core is stabilized by nuclear forces and neutron degeneracy pressure. The 
core is divided into two distinct regions: the subsonic, "homologous" inner core and the 
supersonic outer core. The inner core becomes thermalized and stabilized. 

In the accretion phase, material in the outer core accretes onto the inner core after the 
thermal stabilization of the inner core, leading to the increase in mass of the inner core. We 



1 The stability of the c ore depends on the adiaba tic index. If the average adiabatic index is less than 
4/3, the star is unstable (jShapiro fc Teukolskvlll983l ). In the case of very massive stars, the pressure is 
dominated by radiation pressure, whose adiabatic index is 4/3. Since there are additional components such 
as the non-relativistic nucleon, whose adiabatic index is 5/3, the mean adiabatic index is between 4/3 and 
5/3 so that the core is mildly stabilized. 
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here note that the inner core continues to shrink after stabilization due to the accretion from 
the outer core. These features are seen in Figured] In the top panels and bottom left panel, 
the accretion shock can be seen at ~ 2 x 10 7 cm for t = 2.0 ms and ~ 1.5 x 10 7 cm for 
ms before BH formation. The bottom right panel shows the increasing mass of the inner 
core. Because the outer core is out of sonic contact with its inner counterpart, an accretion 
shock is formed between the two regions 0. It should be noted that the region just behind 
the accretion shock, that is the surface of the inner core, is the main source of neutrinos. 
This can be seen in Figure [2] which shows snapshots of the cooling rate of neutrinos, Q u , 
as a function of the mass coordinate. Comparing with Figure [I] we find that the maximum 
cooling rate is just behind the accretion shock. 

In the case of very massive stars, the inner core is more massive than the maximum 
mass of neutron stars so that these stars collapse to black holes directly. As for the non- 
rotating model, the initial mass of the BH is ~ 5M . This is signif i cantly smaller than 
our previous result, 20M Q for the non-rotating model in ISuwa et al.l (j2007bl ). because of 
general relativistic effects, newly implemented in this study. Although the core is stable in 
Newtonian gravity, the hot core is not stable when general relativistic effects are considered 
because the thermal energy also contributes to the gravitational potential. In our case the 
temperature is so high that general relativistic effects become important because the critical 
temperature, a bove which the core become s unstable, is T cr ; t ~ 1O 12 K(1OM /M), where M is 
the core mass (IShapiro &; Teukolskylll983l ). Our results show the inner core mass is ~ 10M Q 
and the central temperature is ~ 10 12 K (see the discontinuity of bottom right panel in 
Figured] for core mass and top right panel for temperature). This instability drives further 
collapse of the inner part of the inner core, which can be seen in t he radial velocity profil e 
(see again bottom panels in Figured]). This feature is also seen in iNakazato et al.l (120061 ). 
which is a fully general relativistic simulation (see Figure 2 in their result). Although the BH 
is smaller than the inn er core at fir s t, the i nner core wi l l be soo n swallowed completely. Here 
we note that following iFryer et al.l (120011 ); ISuwa et al.l (l2007bl ). BH formation is ascribed to 
the condition 6G ™( r ) > r, where c, G, m(r) are the speed of light, the gravitational constant, 
and the mass coordinate, respectively. After formation, we excise the region inside and then 
treat it as an absorbing boundary. Then, we enlarge the boundary of the excited region to 
take into account the growth from the mass infalling into the central region. Although it is 
not accurate at all to refer to the central region as a BH, we adhere to this simplification in 
order to follow and explore the dynamics later on. 



2 This accretion shock is similar to the stalled and hydrostatic shock front in the framework of ordinary 
core collapse supernovae. The difference is that there is no outgoing shock wave for very massive stars due 
to the softness of the equation of state inside the inner core. 
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Next, we mention the neutrino emission and its relation to dynamics. The time evolution 
of the neutrino luminosity, L v = J Q u dV, is shown in Figure [31 As the collapse proceeds, the 
density and temperature of the central both increase, increasing the neutrino luminosity. As 
the density increases during core collapse, the opacity rises and the neutrinosphere is formed, 
in which neutrino emission-absorption equilibrium is achieved. The relation of the radii of 
the inner core surface, R 1C , and neutrinosphere, _R„, is very important for the luminosity of 
neutrinos because the inner core is the main source of neutrino emission as already shown 
in Figure [2J Roughly speaking, if R ic > R u , the emissivity is the largest just behind the 
shock so that the luminosity is also large. On the other hand, the luminosity becomes 
small if R v > R\ c because the emitted neutrinos can not escape directly and interact with 
matter. In this case, the maximum value of the emissivity appears between the inner core 
surface and the neutrinosphere. In Figure H] we show the evolution of radii of the inner 
core, neutrinospheres, and black hole as a function of time. The inner core surface goes 
down inside the neutrinospheres, which results in the peak in the evolution of luminosity 
shown in Figure [3] (see dots A and B in Figure [3] and [4]) . In addition, after BH formation 
the luminosity decreases due to the absorption of neutrino emitting matter. Finally, the 
neutrinosphere gets swallowed by the BH and the emission of neutrinos almost terminates. 



4. Rotational Collapse 

4.1. Features of Neutrino Luminosity 

We now investigate how rotation affects the properties of the neutrino emission. Peak 
luminosity and total energy emitted by neutrinos for different values of the initial rotation 
rate, T/|W|init, are shown in Tabled] Interestingly and maybe against the intuition, it is 
found that the peak neutrino luminosity increases with the initial rotation rate (see Table [1]). 
We point out that this is because the inner core can partly extend outside the neutrinosphere 
due to rotation so that R ic > R v is achieved, and the degree of the deformation becomes 
larger for more rapidly rotating models. In the following, we explain this feature in more 
detail. 

Figures H] and [5] show the time evolution of radii of neutrinospheres and inner core for 
non-rotating model and rotating model (T/|W|hut = 0.5%), respectively. Note again that 
we set the model with T /|W|init = 0.5% to be the reference rotating model according to 
Heger et al.l (120001 . 120031 ). For the non-rotating model (Figure HJ), the inner core becomes 
smaller than the surfaces of all of the neutrinospheres just before BH formation (see evo- 
lutions after B in Figur^. On the other hand, for the rotating model, the inner core is 
broadened in the equatorial plane due to the centrifugal forces (see left panel of Figure [5]) 
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and exists outside of the neutrinosphere of vx- Moreover, the inner core surface is located 
outside of the neutrinosphere of v e along the polar axis (right panel). The above effects make 
the neutrino luminosities larger for the rotating models. In Figure El the global shapes of 
neutrinospheres of v e and Ux, and the inner cores for spherical model and the rotating models 
with r/|W|init = 0.1%, and 0.5% are presented. It is clearly seen from the comparison with 
the non-rotating model that neutrinospheres are more deformed for more rapidly rotating 
models. As the initial rotation rate becomes larger, the deformed inner core is shown to 
have larger surfaces approaching the neutrinospheres than the non-rotating model. As for 
the model with T/|iy|i nit = 0.5%, the inner core extends outside the neutrinosphere of vx 
entirely. 

The time evolution of the neutrino luminosity for models with different initial strength of 
rotation are shown in Figure [7J For slow rotation models (T/|iy| init < 0.1%), the luminosity 
properties of v e and v e are shown not to depend on the strength of rotation. On the other 
hand, the luminosity of vx gets larger as the initial rotation gets faster. This is because 
the radius of the neutrinosphere of vx is smaller than those of v e and v e so that vx is most 
sensitive to the rotation. The inner core surface, which is in the neutrinospheres for the 
non-rotating model, is enlarged by rotation and falls close to the neutrinosphere of vx on 
the equatorial plane, which leads to the increase of Vx luminosity. In fast rotation models 
(Ty|W|hut ^ 0.3%), the luminosity of v e and v e also get larger with rotation for the same 
reasons as for the vx case. The luminos ity of vy, however, gets smaller as the rotation gets 



faster. This feature was already seen in lSuwa et al.l (l2007bl ). This is because rapid rotation 



suppress the contraction of the most inner part causing the density and temperature not to 
increase as much as in the case of slower rotation. Therefore, the emissivity of neutrinos also 
does not increase very much. This suppression affects the v e and v e emission less because 
neutrinospheres of these species are located at larger radii. 

From Figure [7J we can see drastic decreases in neutrino luminosities near the epoch of 
black hole formation, particularly for the rapidly rotating models. The slope of the neutrino 
luminosity depends on the relative position of the neutrinospheres and the inner core. For 
slow rotation models, only vx shows such a decrement because the neutrinosphere of vx is 
the smallest so that it is absorbed first and is inside of the BH when the BH forms. On 
the other hand, the neutrinospheres of v e and v e remain outside of the BH at the onset of 
BH formation and are absorbed as the BH grows. The luminosity of v e and v e decrease 
with the dynamical timescale (~ O(10) ms) of the region of the neutrinospheres. For fast 
rotation models, all species show the drastic decrease because all neutrinospheres are partly 
or entirely in the BH region when it forms (see Figure Q again ) because the rapid rotation 



leads to a larger initial mass of the BH (see lSuwa et al.ll2007bl ) 
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We find that rotation enhances the total energy released by neutrinos. The total energy 
emitted in the form of neutrinos during the collapse is summarized in Table [TJ It can be 
seen that for the non-rotating model the total radiate d energy amounts to 6 .6 x 10 53 ergs. 
This result is in good agreement with the simulation of iNakazato et al.1 (120061 ). who included 
neutrino transport during the collapse of a 300 M Pop III star. They found a total energy 
output of ~ 4 x 10 53 ergs. The stronger rotation models show larger energy emission. This is 
because the rotation enhances the neutrino luminosity as already described. In addition the 
rotation delays BH formation by the centrifugal force, which leads to a longer duration of 
neutrino emission. These two effects account for the increase of radiated energy by neutrinos 
as the strength of rotation increases. 



4.2. Features of Neutrino Energy 

We discuss the effects of rotation on the emitted neutrino energies in this subsection. 
We consider the averaged neutrino energy, (E u ) = f Q u dVdt/ f N u dV, where N u is number 
density of emitted neutrinos. Figure [8] shows the difference of the average energy as a 
function of T/|W|mif It can be seen that the average energies increase with the initial 
rotation rate for T/|W|init < 0.1% irrespective of the flavors, which is basically the same 
reason as discussed in the previous subsection. It should be noted that the energy of v x is 
more sensitive to the rotation strength than the other species because the neutrinosphere 
of v x forms the deepest in the core. In fact, only v x is affected for relatively slow rotation 
(T/|V^|i n it < 0.1%). As the rotation becomes stronger, the energies of v e and v e are gradually 
changed. Although the average energies of neutrinos b asically foll o w the standard order 



sequence, namely (E Ue ) < (E 9e ) < (E Ux ) as described in llocco et all (120051 ). (E Ux ) is found 
to be smaller than (E Pe ) for the model with T/|lF| init = 2%. The strong centrifugal force 
interrupts the contraction of the inner part of the core, leading to the suppression of the 
releasable gravitational energy and thus smaller energy of neutrinos. 

Next, we focus on the anisotropy of the neutrino energies and emission. The spatial 
distribution of the average energy of v e just before BH formation, with and without rota- 
tion are depicted in Figure [9j It can be seen that due to the rotational flattening of the 
neutrinosphere (blue line) near the polar regions (bottom panel), the neutrinos achieve the 
higher temperature just inside the neutrinosphere, leading to the emission of higher energy 
neutrinos. In addition, the direction of anisotropy of the energies can be seen in the rotating 
model. In the polar direction, the inner core surface is almost coincident with the neutri- 
nosphere (see Figure [5] and O again) so that neutrinos with high energy produced just inside 
the inner core can escape. On the other hand, the optical depth is ~ 5 at the inner core on 
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the equatorial plane and few neutrino with high energy can escape. These effects, irrespec- 
tive of the neutrino flavors, make the neutrino spectrum of the rotating model broader and 
harder than the non-rotating model, which will be useful for the discussion in the following 
subsection. 



4.3. Neutrino Spectrum and Relic Background 



In this subsection, we discuss the features of the neutrino spectrum and the effects of 
rotation. We calculate the neutrino luminosity spectrum as follows: 



dL,_ 



de> 



-(*) 



dVQ Vl 



dP u 

~dV 



(3) 



where Q Ua is the emissivity for u a ( 
Fermi-Dirac distribution, 

dP v 



u e ,u e , and vx) and dP Va /de is the normalized 

(4) 



J2 



7T 4 is 



de' 3C 3 T3e^ + r 

where e' is the energy of emitted neutrinos in the source frame and T a = 180^3 (t' Ua ) /7 
the effective neutrino temperature with \^ Vt j being the local average energy of v a . Figure [TU1 
shows the time integrated energy spectrum, J dt(dL Ua /de'), of u e , v e and vx for the model 
with T/lWlinit = 0.5%. It can be seen that the spectrum of vx (dashed line) is harder 
than v e (solid line) and v e (dotted line) as already described in the previous subsection. 
Figure [TT1 shows time integrated energy spectrum of u e for models with 2^/1 W^|i n it = 0,0.5 
and 2%. It is obvious that the spectra of rotating models (dashed and thick dotted lines) 
are higher and harder than non-rotating model (solid line). In this figure, the thin dotted 
line represents the single energy spectrum. We use the total energy emitted by z/ e , J L Ue dt, 
and total average energy, (E Ue ), of the model with T/|VT|i n i t = 0.5% for this spectrum 
to compare the spectrum of the model with the same rotation strength. By making the 
comparison between thin and thick dotted lines, it can be seen that the calculated spectrum 
(thick dotted line) is harder than one with a single energy. This is because the neutrino 
temperature is not a simple, single temperature component, but rather consists of multiple 
temperature components, depending on the time and position. It should be noted that 
the time evolution of the average energy from Pop III stars is drastic unlike the ordinary 
core-collapse supernovae, which do not show significant change of the average energy on the 
cooling phase of proto-neutron stars (see, e.g.. lLiebendorfer et al.ll2005l ). so that the spectrum 
becomes broade r than the ordinary supernovae. Th e single energy spectrum employed in 
previous works (jlocco et all 120051 ; iDaigne et al.l 120051 ) should have led to underestimates in 
the neutrino number in the high energy regime. 
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We are now in a position to discuss the relic neutrino background from Population III 
stars. For si mplicity, we assum e that all Pop III have identical emission characteristics. 



According to llocco et al.l (120051 ) the differential flux expected at Earth can be written as 



cfnm^^L r dz (l + z) ^ z )^ ( 5 ) 



de J 11 M IH Jo de 

where e = e'/(l + z) is the observed energy of neutrinos, fm is the fraction of all baryons 
going through Pop III, n 7 ~ 410cm -3 is the cosmic microwave background photon density at 
redshift zero, rj ~ 6.3 x 10~ 10 is the cosmic baryon-photon ratio, mjy is the nucleon mass, Mm 
is the mass of Pop III stars, ip(z) is normalized star formation rate, and dN Ua /de' represents 
the number spectra at the source frame, which is also calculated by volume integration with 
the normalized Fermi-Dirac distribution (Eq. (j4])). We assume that the star formation 
rate is concentrated around a redshift z ~ 10, that is, ip{ z ) = $( z ~ 10)- I n addition, we 
assume that all Pop III stars have the same mass, M in = 3OOM . We used the result of 
the T/|W|i n it = 2% model to calculate the relic neutrino flux. Since star formation history 
is not well constrained observationally at high z, we consider samples of //// as 10 _1 , 10~ 3 , 
and 10" 5 . Fi gure [12] shows the calculated number flux spectrum of v e . If we adopt the 
largest value for fjjj, the diffuse anti-electron neutrinos from Pop III dominates those from 
ordinary supernovae (dashed line) below ~ 7 MeV. The dotted line represents the flux from 
the model with T/|iy|j n it = 0% and fm = 10~ 3 . Comparing with the second solid line from 
the top, we can see the amplification of the relic neutrino flux by the stellar rotation. 



5. Summary and Discussion 

We performed a series of 2D hydrodynamic simulations of rotational collapse of Pop 
III stars. We computed 8 models, changing the initial strength of the angular momentum 
in a parametric manner. We studied systematically the deformation of neutrinospheres and 
investigated the degree of anisotropic neutrino emission from the point of view of the neutrino 
luminosity and the average energy, and their impacts on the neutrino spectrum. Then we 
found the following, 

1. The neutrino luminosity increases for more rapidly rotating cores. This is because the 
inner core, which is the dominant neutrino source, is deformed due to the rotation 
and can extend outside the neutrinospheres. In particular, vx is most sensitive to the 
rotation because its neutrinosphere is located at the smallest radius and nearest to the 
inner core. The total energy emitted by neutrinos also increases with the strength of 
rotation. 
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2. The average energy of emitted neutrinos gets larger as rotation becomes faster. In 
addition, the energy is larger in the polar direction than in the equatorial plane due to 
the rotation. The neutrino spectrum gets broader and harder for rotating models. As 
a result, the number spectrum of the relic neutrinos also gets harder so that the relic 
neutrinos from Pop III stars overwhelm those from ordinary core-collapse supernova 
below ~ 7 MeV, provided a fraction //// = 0.1 of all baryons forms Pop III stars. 

As for our numerical computations, we shall address a few of our assumptions. In this 
study, we mimicked the neutrino transfer by the leakage scheme. Although the scheme is 
a radical simplification, we checked that in the case of spherical collapse, the obtained in- 
tegral quantities such as the to tal energy emitt e d by neutrinos and their average energy 
are very close to the results by iNakazato et al.l (120061 ). in which the Boltzmann neutrino 
transport was solved. We ignored the angle dependence of the neutrino emission, assum- 
ing that the neutrino track is radial. This approximation will break down just before BH 
formation because the neutrinosphere becomes oblate and the neutrino emission would be 
highly aspherical. Nonetheless, the multi-angle and multi-dimensiona l neutrino transfer is 
still computationally prohibitive (recently reported by lOtt et al.ll2008l in the framework of 
core-collapse supernova). The angle dependence is expected to make the neutrino spec- 
trum harder. This is because the neutrinos produced in the vicinity of the equatorial plane 
might be able to escape to the polar direction even if they cannot escape radially, thus it 
should be taken into account. Furthermore the Newtonian simulation in this paper, albeit 
with the general relativistic corrections, is nothing but an idealized study to describe the 
dynamics with BH formation. It is by no means definitive. In order to expl ore these phe 



nomena in more detail will require fully g eneral relativistic simulations (e.g.. IShibatal 12003 



Sekiguchi fc Shibatal 120051 ; iLiu et al.l 120071 . for BH formation of massive stars) but with the 
microphysical ingredients at least at the same level of this study, which is beyond the scope 
of this paper. 

Bearing these caveats in mind, we state some discussions and speculations about the 
relic neutrinos from Pop III stars. In estimating the relic neutrino background, we employed 
a parameter fuj, the baryon fraction of Pop III stars. Because this parameter is entirely 
unknown so far, we employed three values of 10 _1 , 10 -3 , and 10~ 5 . The largest value is 
obtained for the explanation of the infrared b ackground excess by UV photons from Pop 
III stars ( jSantos et al.l 120021 ; iDwek et al.l 120051 ). So our estimation should be interpreted 
as giving an upper bound. As for the star formation history, we employed only the delta 
function one because the baryon fraction dominates the uncertainty a nd the other ones 
such as the evolution of star formation rate are minor corrections (see, Iocco et al. 2005 



Daigne et al.l 120051 ; INakazato et al.l 120061 . for the discussion of the dependence on the star 
formation history). In addition, we assume that all Pop III stars have the same mass of 
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300M Q for simplicity. The dependence of the progenitor masses should be also clarified, 
although the qualitative features induced by rotation studied here will not change. 

The predicted foregrounds for the relic neutrinos in the relevant energy range computed 
here are the neutrinos produced in the Sun and the anti-neutrinos from nuclear reactors. 
These ( ant i-) neutrinos, however, could be removable in principle because the distribution 
of these foregrounds are not isotropic. Other possible foregrounds, w hich are isotropic, ar e 



extra-galactic neutrinos produced by thermonuclear reactions in stars (jPorciani et al.ll2004l ). 
On the other hand, as for anti-neutrinos the contribution of Pop III stars are dominant in 
the range below a few MeV. Thus we speculate that the detection of this signature could be 
the first direct probe of this as yet unseen metal-free first generation of stars. 
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Fig. 1. — Time evolution of density (top left), temperature (top right), and radial velocity 
(bottom left) for non-rotating model as a function of radius in 10 7 cm. In addition, the radial 
velocity as a function of mass coordinate in M & (bottom panels) is presented for comparison. 
The thin-solid, dotted, dashed, dot-dashed, and thick solid lines are snapshots at 24.1, 10.1, 
2.0, 0.9, and ms before black hole formation, respectively. 
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Fig. 2. — Snapshots of total neutrino cooling rate Q v in 10 32 erg s _1 cm -3 as a function of 
mass coordinate in M , for the non-rotating model, showing that the most efficient cooling 
appears just behind the accretion shock. Times are relative to black hole formation time. 
Compare this figure with Figure [I], which shows the hydrodynamical evolution for this same 
model. 
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Time after black hole formation (s) 



Fig. 3. — Time evolution of neutrino luminosity L v in 10 erg s for the non-rotating 
model. Solid, dotted and dashed lines represent v e , v e and u x , respectively. Circles A and 
B correspond to the peak luminosity for v e and vx- 
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Fig. 4. — Radii of neutrinospheres of v e (solid), v e (short-dashed), ux (dotted), inner core 
surface (long-dashed) and BH (dot-dashed) versus time after black hole formation in seconds 
for the non-rotating model. The radii of the v e and v e neutrinospheres are similar because 
the electron fraction does not reach small value as in the case of core-collapse supernova. The 
inner core is inside all of the neutrinospheres just before BH formation so that the luminosity 
of the non-rotating model does not get large in Figure [7J 



-21- 



(a) Equatrial Plane 
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(b) Polar Axis 
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Fig. 5. — Time evolution of radii for a rotating model with T/|ty| init = 0.5%. The left panel 
is in the equatorial plane and the right panel is along the polar axis. The line types are the 
same with Figure H] except that the neutrinosphere of v e is not shown in this figure because 
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Fig. 6. — Positions of neutrinospheres and the inner core. The first panel shows the non- 
rotating model, the second T/|H^| init = 0.1% and the third T/|W|init = 0.5%. The solid, 
dotted and dashed lines give the locus of the neutrinospheres of u e , vx and the inner cores, 
respectively. 
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Fig. 7. — Time evolution of neutrino luminosity for all models. v e , u e and vx are represented 
by solid, dotted and dashed lines, respectively. 
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Fig. 8. — The total average energy of neutrinos, (E u ) = f Q u dVdt/ f N u dV, where N u is 
the number density of emitted neutrinos, as a function of T/lVKlmit- The solid, dotted and 
dashed lines represent u e , v e and ux, respectively. 
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Fig. 9. — Average energy distribution of v e and the shape of neutrinosphere just before 
BH formation for models with T/|jy|i nit = 0% (left) and 0.5% (right). The blue solid 
line represents the position of the neutrinosphere. The central black part corresponds to 
the region with the optical depth of neutrino being beyond 10, where neutrinos are almost 
trapped. 



Table 1: Results of the numerical simulations. 





r peak 

V 


-E'total 


(%) 


(10 54 erg/s) 


(10 53 er 





1.11 


6.57 


0.01 


1.17 


6.78 


0.05 


2.17 


8.30 


0.1 


4.19 


11.3 


0.3 


7.55 


22.2 


0.5 


8.22 


29.6 


1 


8.64 


40.8 


2 


8.44 


52.4 



Note. — T/|W|init means the initial ratio of the rotational energy, T, and the gravitational energy, W . 
jjpeak represents the peak total luminosity of neutrinos. -Etotai is the total energy emitted by neutrinos. 
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Fig. 10. — Energy spectra of neutrinos for the model with T/lVKlinit — 0.5%. The red solid, 
green dotted and blue dashed lines represent u e , v e and ux, respectively. 
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Fig. 11. — Energy spectrum of v e for model with different values of T/|W|init- The red solid, 
green thick-dotted and blue dashed lines represent T/|W^| ini t = 0, 0.5 and 2%, respectively. 
The green thin-dotted line shows the spectrum with the single temperature Fermi-Dirac 
distribution function (see §4.31 for discussion). 
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Fig. 12. — Calculated number flux of v e (red solid lines). The employed baryon fraction of 
Pop III stars, fin, is 10 _1 , 10 -3 and 10 -5 from top to bottom. As for the solid lines, we 
employed the result of the T/|W|init = 2% model. The dotted line represents the flux with 
the T/|W| init = 0% model and //// = 10~ 3 . For comparison we show the diffuse flu x from 
ordinary core-collapse supernovae (blue dashed line) according to lAndo et al.l (120031 ). 



